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The evolution of resistivity versus temperature (p(T)) curve for different doping elements, and in 
. the presence of various defects and clustering are explained for both diluted magnetic semiconduc- 

• tors (DMS) and manganites. Here, we provide unambiguous evidence that the concept of ionization 

' energy (Ei), which is explicitly associated with the atomic energy levels, can be related quantita- 

/"S , tively to transport measurements. The proposed ionization energy model is used to understand how 

the valence states of ions affect the evolution of p(T) curves for different doping elements. We also 
0^ ' explain how the p(T) curves evolve in the presence of, and in the absence of defects and clustering. 

The model also complements the results obtained from first-principles calculations. 
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C3 . I. INTRODUCTION 

(— > 

Ferromagnets have the tremendous potential for the development of spintronics and subsequently will lay the 
foundation to realize quantum computing. The field of spintronics require the incorporation of the spin-property 
of the electrons into the existing charge transport devices Parallel to this, the technological potential of DMS 
(Rcf. 0) is associated to spintronics-device development, whereas manganites that show a large drop of resistance 
below Tc lead to the colossal magnetoresistance effect (CMR), which is also important in the new technologies such 
Q as read/write heads for high-capacity magnetic storage and spintronics As such, applications involving both 
DMS and manganites very much depend on our understanding of their transport properties at various doping levels 
and temperatures (T). In addition, DMS also has several interesting physical properties namely, anomalous Hall- 
effect [|[, large magnetoresistance in low dimensional geometries Q, the changes of electron-phase-coherencc time 
in the presence of magnetic impurities Q and negative bend resistance Q- As for the transport properties, there 
are several models developed to characterize the resistivity of DMS. In particular, the impurity band model coupled 
with the multiple exchange interactions for T < Tc for Gai-aMn^As was proposed @- The electronic states of the 
■ impurity band can be either localized or delocalized, depending on doping concentration or the Fermi-level {Ep). If 
Ep is in the localized-state, then the conduction is due to carrier hopping. If Ep is in the extended-state, then the 
conduction is metallic and finite even for T = (Ref. 0). On the other hand, the spin disorder scattering resistivity as 
f-*; | a function of magnetic susceptibility can be used to estimate the magnitude of J ex (the ferromagnetic (FM) exchange 
'^""-i , interaction energy) Q. Moreover, there are also theories that qualitatively explain the conductivity for T > 0, namely, 
the Kohn-Luttinger kinetic exchange model and the semiclassical Boltzmann model [ll|, EH 13 • 

Apart from that, for manganites, the one- and two-orbital models fl4j and the phase separated resistivity model [Tol 
have been used to qualitatively describe the resistivity curves for T > 0. However, in all these approaches, we 
are faced with two crucial problems, the need (i) to explain how the resistivity evolve with different doping elements, 
without any a priori assumption on carrier density and (ii) to understand how defects and clustering affect the 
evolution of p(T) curves. Here, we show unequivocally, a new method to analyse the evolution of p(T) curves 
for different doping elements using the concept of the Ei invoked in the Hamiltonian and Fermi-Dirac statistics. 
In doing so, we can also understand the evolution of p(T) curves in the presence of defects and clustering, which is 
important for characterization of spintronics devices. The Ej concept has broad applications, where it has been applied 
5_i ' successfully for the normal state (above critical temperature) of high temperature superconductors [l?], EH, Oil an d 
ferroelectrics [20j]. The Ej model is for compounds obtained via substitutional doping, not necessarily homogeneous 
or defect-free. 
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II. IONIZATION ENERGY MODEL 
A. Carrier density 

A typical solid contains 10 23 strongly interacting particles. Therefore, their universal collective behavior is of 
paramount interest as compared to the microscopic details of each particular particle and the potential that surrounds 
it. This universal collective behavior, being the focal point in this work, arises out of Anderson's arguments in More 
is Different. [2l[ That is, we intend to justify a universal physical parameter that could be used to describe the 
association between the transport-measurement data and the fundamental properties of an atom. In view of this, we 
report here the existence of such a parameter through the Hamiltonian as given below (Eq. ([T])). The parameter is 
the ionization energy, a macroscopic, many-electron atomic parameter. 

^=(^±0^ (!) 

where H is the usual Hamilton operator and Eq is the total energy at T = 0. The + sign of ±£ is for the electron 
(0 — > +oo) while the — sign is for the hole (— oo — > 0). Here, we define the ionization energy in a crystal, £ = SJ eal is 
approximately proportional to Ei of an isolated atom or ion. Now, to prove the validity of Eq. {1} is quite easy because 
£ is also an eigenvalue and we did not touch the Hamilton operator. Hence, we are not required to solve Eq. ([I]) in 
order to prove its validity. We can prove by means of constructive (existence) and/or direct proofs, by choosing a 
particular form of wavefunction with known solution (harmonic oscillator, Dirac-delta and Coulomb potentials) and 
then calculate the total energy by comparison. In doing so, we will find that the total energy is always given by 
Eq ± £, as it should be (see Appendix and Ref. [l^)- For an isolated atom, the concept of ionization energy implies 
that (from Eq. 

±£ = Skin — E + Vpot = ±Ei, 

where Ei is the ionization energy of an isolated atom. The corresponding total energy is Eq ± £ = Skin + Vpot = 
EqzszEj. Whereas for an atom in a crystal, the same concept of ionization energy implies that ±£ = Skin — Eq + V po t + 
Vmany-body = Ej + Kiany-body = ±SJ cal . Here, V^any-body is the many body potential averaged from the periodic 
potential of the lattice. The corresponding total energy is Eq ±£ = Ski n + Vpot + V man y-body = Eq ±Ej + V man y-body = 
Eq ± Ej ca[ . Here, Sj cal is the ionization energy of an atom in a crystal. The exact values of Ej are known for an 
isolated atom. That is, one can still use Sj obtained from isolated atoms for r — ► oo in order to predict the evolution of 
resistivity versus temperature curves for different doping elements. Therefore, Eq. ([1]) can be approximately rewritten 
as Hip cx (EQ±Ei)ip. It is obvious from the Hamiltonian given in Eq. JT]) that we cannot use it to determine interactions 
responsible for FM and the magnitude of Tq, simply because we have suppressed the Vmany-body On the other hand, 
if we have a free-electron system, then the total energy equation is given by E — Eq ± YTi ^iTj = So ± /3 J2j Eij 7 
where, Eq is the total energy of the compound at T = and j is the sum over the constituent elements in a particular 
compound. We also define here, (3=1+ , where (V(x)) is the averaged many-body potential value. Here, E\ cai 

can be noted as one of the many-body response functions with respect to transport properties. The total energy can 
be rewritten as 



E =^±EE^ 

* 3 

= Eq ± [Skin — Eq + V po t + Vmany-body] 

= Skin + Vpot + V ro any-body & f° r electrons ± -> + 

= Skin + Vtotal 

= Skin ^ implies free electrons. 

Note that we have substituted for SJ cal since the concept of ionization energy is irrelevant here simply because the 
electrons in these metals are free and do not require excitations from its parent atom to conduct electricity. As such, 
the carrier density is constant and independent of temperature. Whereas, the scattering rate is the one that determines 
the resistivity with respect to temperature, impurities, defects, electron-electron and electron-phonon interactions. 
Therefore, the total energy from Eq. |T|) carries the fingerprint of each constituent atom in a compound and it refers 
to the difference in the energy levels of each atom rather than the absolute values of each energy level in each atom. 
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Hence, the kinetic energy of each electron from each atom will be captured by the total energy and preserves the 
atomic level electronic-fingerprint in the compound. Parallel to Eq. ([T]), the electron and hole distribution functions 
can be derived as [13, EH (see Fig. [T]) 



fe{Ea,Ei) = 



fh(E ,Ei) 



1 

e [{E +E,)-E^\/k B T + i 
1 

e [£< r 0) -(B --E / )]/fc B T + i 



(2) 



where Ep^ is the Fermi level at T = 0, independent of doping concentration and ks is the Boltzmann constant. 
From Eq. @ , one can surmise that large Ei corresponds to the difficulty of the electron to conduct due to a large 
Coulomb attraction between the electron and the ionic core. This effect will give rise to low conductivity. The 
variation of Ej with doping in our model somewhat resembles the variation produced from the impurity band model 
with doping [H, [23| . The difference is that our model fixes the Fermi level as a constant and we let Ej capture all 
the changes due to doping, consistent with ionization energy based Fermi-Dirac statistics (iFDS). Whereas, in the 
impurity-band approach, both the impurity band and the Fermi level evolve simultaneously with doping, consistent 
with Fermi-Dirac statistics. 



FDS, iFDS (T1 > T) 
FDS, iFDS (T = 0) 
FDS, iFDS (T > 0) 




E, (E + Ej 



FIG. 1: Standard Fermi Dirac (FDS) and ionization energy based Fermi-Dirac (iFDS) distributions for temperatures, T 
T > and Tl > T. 



The carrier density can be calculated from 



= / f e {E 0l E!)N e (E )dE a , 



(3) 



where N e (Eo) is the density of states and Ep is defined from now on as the Fermi level at T = so as to comply 
with Eq. Q. 



B. Spin-orbit coupling 

Here, we will show the relationship of the ionization energy to the energy level splitting and the spin-orbit coupling, 
as well as their association with resistivity. The latter association is crucial because spin-orbit coupling is found to 
be an an important phenomenon that influences the electronic properties of the magnetic semiconductors (24l . [25| . 
ft is well established that the energy associated with both the relativistic correction and the spin-orbit coupling for 
Hydrogen-like atoms is given by the fine structure formula, Ef s (Ref. [IH) 
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(4) 



where ctf s is the fine structure constant, k = j + 1/2, j = I ± 1/2, where j is the total angular momentum, / is 
the orbital angular momentum and E n is the non-relativistic energy level. Z and n are the atomic and the principal 
quantum numbers, respectively. Using Eqs. |T]) and ((¥]), we can show the extreme magnitude of the energy level 
splitting (i.e., between k — 1 and k — n) is 
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AE fs = (a fs Z) 2 (E ±On-^^. (5) 

Thus, the magnitude of the energy level splitting (AEf s ) is proportional to the ionization energy (£) and Z 2 , 
while it is inversely proportional to n. That is, AEf s also decreases with decreasing ionization energy (see Fig. [2|. 
Therefore, we can surmise that a system with large ionization energy gives rise to large spin-orbit coupling, which 
is important for spin-injection [27| in p-type DMS. As discussed earlier however, a large ionization energy also leads 
to a small carrier density (from Eq. ©), which in turn implies that spin-orbit coupling competes with the electrical 
conductivity. In addition, Fig. [5] points out that A-E/ S i_>2 and £1^2 are larger than AEf s i—>3 and £2— >3, respectively, 
in accordance with Eqs. (U) and (J5]). Therefore, we can surmise that A_E/ S i_>2 °c £i->2; A-E/gi—^ oc ^2->3, A-E/si—^ 
oc £3^4 and so on. 



£.= 1 Jl 



^ 

U 3 VL 



#l-> 3 



/&!-> 2 



i-» 2 



£=3,7 = 5/2 
k = 2,j = 3/2 
*= l,y = l/2 

k = 2J = 3/2 
£=1,7 = 1/2 



(£„ + On=l=£o 



£ = l,j= 1/2 

FIG. 2: Energy levels of hydrogen-like atom and the energy level splitting (not to scale). £n=i,2,3 is the standard energy level 
notation while, (Eq ± £) n =i,2,3 is the new notation introduced from Eq. ([5}. 



C. Resistivity model based on ionization energy 



In order to derive the resistivity model as a function of ionization energy, we need an expression that connects the 
carrier density with the total current. As a consequence, we propose that the total current in ferromagnets consists of 
contributions from both non-FM and FM phases, which is J = ^2 J v . with v = e^, se^ , sde- 1 , h^. sh^ , sdhK Here, se^ 
and sh^ are the spin-assisted electron and hole, respectively below Tc, influenced by the T-dependent magnetization 
function (M(T,Mq)) through the spin disorder scattering rate (tsd)- ^ and are the electron and hole with 
electron-electron scattering rate (r e _ e )- sde^ and sdh^ are the electron and hole with spin disorder scattering and 
M(T, Mo) = constant (Const.) for T > Tc- For convenience, the spin-up, f denotes the direction of the magnetic field 
or a particular direction below Tc, while the spin-down, J, represents any other direction. That is, a non-FM phase 
gives rise to J Eyh and J sde ,sdh- Whereas, J se , s h originates from a FM phase where (J s de..sdh,Jr..h) -> J S e,sh for T < T c . 
Hence, this resistivity model is fundamentally different from the phase separation model [151 ]. That is, J e ,/i(oc r e _ e ) 
< J seiSfi (oc t S d[M(T,M )]) for T < T c , while J e , h (oL r e _ e ) and J s d e ,sd/i(oc Ts£>[Const.]) contribute for T > T c due 
to the non-existence of the FM phase. In other words, for T > Tc, the whole system is a non-FM phase and its 
conductivity is determined by J e ,/i(oc r e _ e ) and J s de,sdh(oc tsd [Const.]). For T < Tc, some of the non-FM phase 
becomes FM due to spin and its conductivity is determined by J e< h oc r e _ e and J se . s /j(cx tsd[M(T, Mq)]). Simply 
put, we assume that r e _ e and tsd [Const.] contribute in the non-FM phase above Tc- Below Tc, we have both FM 
and non-FM phases contributing to the resistivity via r e _ e and tsd[M(T, Mo)]. As such, the total current can be 

simplified as J = j\ + [Jj e ,Jg de ] = J e + [Jse,Jsde] if the considered system is n-type, while J = Jh + [J s h,Jsdh] if it is 
p-type. J e and Jh are the spin independent charge current due to electron and hole, respectively, and are influenced 
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by T e _ e . J sc ie and J sc ih are also the spin independent charge current due to electron and hole, respectively, but are 
influenced by rg^ [Const.] in the non-FM phase, which is valid only above Tq. J se and J s h are the spin-assisted 
charge current in the FM phase only, influenced by the tsd[M(T, Mq)]. Thus, the total resistivity (for n or p-type) 
can be written as (after making use of the elementary resistivity equation, p = m/ne 2 T) 



P 1 — Pe.h + [Psde,sdh> Pse,s 



*e,h 



(n,p)e 2 T e _ t 
1 



{n,p)e 2 T SD [Const., M(T, M )] 



m *e,h 


-1 

+ 


m *e,h 


_(n,p)e 2 T e - e _ 




(n,p)e 2 TSD_ 



1 -1 



(6) 

where m* h denotes the effective mass of the electron or hole, e is the charge of an electron. The carrier density for 
the electron and hole (n,p) based on iFDS are given by (after substituting Eq. (J2J) into Eq. ([3])) 



k B T 



2irh 2 



3/2 



{mlf 2 exp 



Ep — Ej 



(7) 



p=2 



k B T 
2ttH 2 



3/2 



(m£) 3/2 exp 



-Ep — Ei 



k R T 



(8) 



The spin disorder scattering resistivity as derived by Tinbergen-Dekker is given by [25 



Psd(T < T c ) 



{mt, h f/iN{2E F )V* j2 



7r(n,p)e 2 fi. 4 

2 



S(S + 1)-S 2 (^1 

xtanhf 3r c M TD (D 
\2TS{S+1)M 



S 



M TD (T) 
M 



(9) 



where Mtd(T, Mq) is the Tinbergen-Dekker magnetization function |28j . Mq is the magnetization at zero temper- 
ature. N is the concentration of nearest neighbor ions (for example, the concentration of Mn). h = h/2ir, h denotes 
Plancks constant and S is the spin quantum number. Equation ^ is equal to the theory developed by Kasuya p9j if 
one replaces the term, tanh [3T c M td (T)/2TS(S+1)Mq] with 1. Substituting l/r e _ e = A eJl T 2 (due to the electron- 
electron interaction and A e ^ is the T independent electron-electron scattering rate constant), together with Eqs. ([9]) 
and (J7J) (or ©) into Eq. ©, then one can arrive at 



Pe,se(T) = 



AB exp [{EjTE F )/k B T] 
AT 3 / 2 [M a (T, Mq)}- 1 + BT- 1 / 2 ' 



(10) 



A = 



A, 



2e 2 (m* eh y/ 2 



2nh 2 



3/2 



B 



' SD 



2ml h N(nE F y/ 2 J 2 x 
e 2 hkT 
N(2E F f/ 2 { m l h f/ 2 



j 2 Mt,Mq). 



Note here that we have invoked the strong correlation through the carrier density, which is a function of the real 
ionization energy. Secondly, to account for different spin polarization below and above Tc , we have n = , 1 /m* = 
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+ \/m*^ and r = T^_ e +T SD {J ex ). Notice that the spin-disorder scattering rate, t sd is a function of the exchange 
energy, J ex . Therefore, by redefining the variables in the elementary resistivity equation in accordance with strongly 
correlated effects, we have now derived a resistivity equation (Eq. flO]) ) suitable for strongly correlated ferromagnets. 
Here, M a (T, Mq) is the T-dependent magnetization function for T < Tc- For T > Tc however, M a (T, Mq) = constant. 
The negative and positive signs in Ei =p Ep are for electrons and holes, respectively. Equation (| 1 0[) points out that, 
Pe,se(T) is semiconducting and it only becomes metallic if Ej ^Ep < T and the FM phase sets in. Below Tc, the FM 
metallic phase increases with lowering T and once it achieves the maximum value (saturation) at a much lower T, 
insulating character may set in as a result of E\ =p Ep > T . We define this temperature T crossover which corresponds to 
the transition from FM metallic to insulating below Tc- Hence, it is clear that our model is fundamentally different 
to the phase separation model as pointed out earlier. The emergence of insulating (T crossovor ) character below Tc is 
an intrinsic property based on the ionization energy model, regardless of defect densities and it was first predicted 
for both high-T c superconductors and ferromagnets fl7j . In a recent analysis using the one- and two-orbital models, 
a similar insulating character below Tc was observed from numerical calculations (l4| . in support of our prediction. 
However, the analysis stops short of explaining why the insulating behavior persists even in the "clean limit" with 
no defects (i.e. pure host material with zero interstitial and vacancy defects). Here we propose that T crossove r or the 
insulating behavior below Tc is intrinsic and is associated with the energy levels through the Ej =p Ep parameter. If 
Ei =F E F > T > Tc, then p e>se (T) cx exp[(£ , z =p E F )/T]. In this limit, a large CMR effect could be observed if the 
metallic-ferromagnetism sets in. The empirical function of the normalized magnetization is defined here as 



Mp {T,M ) = l--^. (11) 

Equation (fTTj) is an empirical function that will be used to extract the magnetization curve from the resistivity 
curve via Eq. (fTOj) . In other words, Eq. pip is used to calculate the magnetization curve, after coupling it with 
Eq. (fTOj) . For example, the magnetization curves associated with Tinbergen-Dekker (TD) (28|, Kasuya (K) [2^| and 
resistivity curves (p) are calculated using 



M TD (T,M ) 



= S(S+1 
x tanh 



s- 



(M TD {T) 



V Mo 
3T C M TD (T) 



2TS(S+1)M 



S 



M TD (T) 
Mn 



(12) 



M K (T,M ) 
= S{S + 1)-S 



V Mo 



S 



Mk(T) 
Mo 



(13) 



and Eq. (fTTj) , respectively. That is, Eqs. (|12[) . (fT3"j) and (fTTj) are separately coupled with Eq. (fTTij) to fit the resistivity 



curves and subsequently obtain 



M TD (T) M K (T) 



and 



MJT) 



M , respectively. Consequently, we can compare and analyze 
ex P (T)/M data. Recall that, a = K (calculated from Eq. (fT3j) ~). 



M ' M 

M (T)/Mq with the experimentally measured, M, 

TD (from Eq. (|12p ). p (from Eq. (flip ), exp (Experimentally determined magnetization curves) 
will drop the term Mq from the equations that contain ^ ' . 

For the electron- or hole-doped strongly correlated non-FM semiconductors, one needs Eq. 
is again based on iFDS, [13, El, [H, H3| 



For simplicity, we 
T4j) given below, which 



Pe,h 



AT 1 ' 2 



cxp 



k B T 



(14) 



where the — sign is for the electron and + sign is for the hole. Equation (fTTj) will be used to justify the importance 
of the term J se even if the resistivity is semiconductor-like in the FM phase. 
The values of Ej can be averaged using the approximation given below 



E 



cal 



. z ' 



(15) 
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where (3 is a function of the many-body potential, Vmany-body and varies with different "background" atoms (host 
lattice). The i = 1, 2,..., z represent the first, second, ...etc ionization energies and, z is the oxidation number of a 
particular ion. Ej cal is actually equal to the energy needed to ionize an atom or ion in a crystal such that the electron 
is excited to a distance r. On the other hand, E\ corresponds to taking that particular electron to r — ► oo with 

^many— body 0- 



o 

C 2 

(1 1 ^ 



□ 



Ol 



□ 



L> ■ A 



K) 



^ ♦Ca 
■ Mn 

OGe 
□ La 

\ mm ^ Ga . 

5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 
Average ionization energy (eV/atom) 



FIG. 3: Average ionization energy distribution for different isolated atoms and valence states. Calculated from Ref. [3jJ using 
Eq. (p)l 

Prior to averaging, the ionization energies for all the elements mentioned above were taken from Ref. [3l| . also 
given in Fig. [3] We can also use Eq. (|T6|) that originates from Eq. (fl~5|) to predict the change in the valence state of a 
particular ion. 



x z+ i i z i 9 

- E ^ + i E^ = 1 E^ 



(16) 



i=z+l 



Here, j Ej=2+i -^fi nas * = z + z + .7 an( l J = 1> 2,.... It is solely due to the multivalencc ion. That is, 
the first term is due to Mn 4+ ion's contribution (Mn 3+ — electron — > Mn 4+ = 51.200 eV atom -1 ), hence j equals 1 
(=4 — 3) in this case and 8 represents the additional contribution from Mn 4+ . The second (i = 1, 2,..., z) and last (i 
= 1, 2,..., q) terms respectively are due to Mn — 3(electrons) — > Mn 3+ and Ga — 3(elcctrons) — > Ga 3+ . Recall that 
q = z = 3+ and i = 1, 2,... represent the first, second,... ionization energies while j = 1, 2,... represent the fourth, 
fifth,... ionization energies. 



III. IONIZATION ENERGY MODEL APPLIED TO FERROMAGNETS 



A. Gai-^MrizAs 

1. p(T>T c ) 

We will first apply the i^j-based carrier density to the resistivity (p{T)) curves of DMS above Tq by means of 
p(T > Tq) oc l/n, in order to avoid the influence of fitting parameters in our analysis. Subsequently, we will discuss 
the evolution of resistivity curves with appropriate fittings. The Gai-^Mn^As system is the most studied material 
among the DMS with Tc around 170 K (Ref. [111 ). T he origin of ferromagnetism and the transport properties in this 
class of materials are still unclear quantitatively [ll| . The evolution of the resistivity curve {p{T,x)) for sample XI, 
Gai-^Mn^As (see Fig. 3 of Ref. d|) is such that p(T, 0.015) > p(T, 0.022) > p(T, 0.035) > p{T, 0.043) > p(T, 0.053). 
That is, the p(T, x) curve shifts downwards with increasing dopant concentration, x. This is expected from the iFDS 
model since Ej for Mn 3+ (18.910 eV atom -1 ) is less than Ej for Ga 3+ (19.070 eV atom -1 ), assuming Mn 3+ substitutes 
Ga 3+ . However, the resistivity curve, p(T, 0.071) is above p(T, 0.035), unexpectedly. 

For sample Y2 (see Fig. 1 of Ref. [32|), p(T, 0.015) > p(T, 0.022) > p(T, 0.035), also complies with the iFDS model. 
Again, unexpectedly we have, p(T, 0.043) < p(T 7 0.053) < p(T, 0.071), where p(T, x) shifts upwards with increasing x. 
That is, p(T,x) switched over from decreasing (expected) to increasing with x at critical concentrations, x c \ = 0.071 
and x C 2 = 0.043 for samples XI and Y2, respectively. These switch-overs that seem to violate iFDS can be explained 
if we can understand what causes E\ to deviate from its averaged value. The only reason that could give rise to such 
deviation is the change in the ions valence states. 



8 



Using Eq. dHJ), we obtain: [<5£/(Mn 4 +)] + i[£/(Mn 3 +) + E I {Mn 2 +) + £?/(Mn 1+ )] = ±[£/(Ga 3 +) + E I {Ga?+) + 
£/(Ga 1+ )] and (51.200 x S) + 18.910 = 19.070, therefore S = 0.003. Since the average Ej (18.910 eV atom- 1 ) 
for Mn 3+ is less than the average Ej (19.070 eV atom -1 ) for Ga 3+ , and if the resistivity curve shifts upward with 
Mn substitution, then z + 5 gives the minimum valence number for Mn (Mn > ' z+i ' + ) which can be calculated from 
Eq. (|16| . If however, the resistivity curve shifts downward with Mn substitution, then z + 5 gives the maximum 
valence number for Mn (Mn < ^ z+S ' >+ ). Consequently, the maximum valence state for Mn in Gai-^Mn^As is 3.003+ 
(Mn <3 003+ ) for the case where the p(T, x) curve shifted downwards with x. If the valence state of Mn is larger than 
3.003+, then p{T,x) is expected to shift upwards with x. Therefore, the switch-over from decreasing to increasing 
p(T,x) with doping, x is due to larger average Mn valence state, and it is calculated to be larger than 3.003+. Next, 
we need to understand what can cause this change in the average Mn valence state. The most likely reason comes 
from Mn occupation at non-substitutional sites, i.e., the Mn 3+ ions do not substitute Ga 3+ ions. This has been 
predicted experimentally where the formation of Mn interstitials (Mn/) is found to be substantial above a critical 
concentration, x c (Ref. Q). Interestingly, occupations at interstitial sites is found to change the charge states of Mn 
ions, as predicted from the first-principles calculations, where Mn 5+ and Mn 4+ are more stable at interstitial sites, 
while Mn 3+ ' 2+ ' 1+ are stable at substitutional sites 33]. In addition, interstitial in the presence of clustering gives rise 
to larger positive charge states as compared to clustering due to substitutional alone in which, larger positive charge 
implies larger valence state for Mn ions [33]. Therefore, we can understand that at certain doping, the average Mn 
valence states increased due to Mn/ and clustering and eventually reduces the carrier density and shifts the p(T, x) 
curve upwards with doping. 



2. p(T > 0) 

Now we will apply the ionization energy incorporated resistivity model to the resistivity measurements [H and the 
fits based on Eqs. (fTP)) and (|14p are shown in Figs.[3|a) and0|b) respectively for Gai-^Mn^As. One needs two fitting 
parameters (A and Ej =F E F ) for p(T > T c ) and another two (B and M p (T,M )) for p(T < T c ). All the fitting 
parameters are listed in Table fl] 



TABLE I: Calculated values of A, B and the ionization energy (£/). "Ann. T" denotes the annealing temperature and T cr 
= T cro3SOVer , is the temperature that corresponds to the transition from FM metallic to insulating character. "Calc." denotes 
calculated values from the fits. 



Samples 


Ann. T (H) 


A 


B 


Ei T E F (Calc.) T C (T cr ) 




°C (Tesla) 


(Calc.) (Calc.) 


Kelvin (meV) 


Kelvin 


Gao.94oMn .o6oAs (a) 


370 (0) 


4.5 


400 


8 (0.69) 


50 (10) 


Ga .93oMn .o7oAs (a) 


As grown (0) 


9.2 


400 


12 (1.04) 


45 (12) 


Ga .93oMno.o7oAs (a) 


370 (0) 


0.02 




280 (24.2) 




Ga .93oMno.o7 As (a) 


390 (0) 


0.03 




400 (34.5) 




La .9Cao.iMn0 3 (ft) 


-(0) 


10 


0.65 


1400 (121) 


222 (-) 


Lao.8Cao.2Mn0 3 (i,) 


- (0) 


10 


1.2 


1300 (112) 


246 (-) 


La .8Cao.2Mn0 3 (f,) 


- (6) 


5 


3.2 


900 (78) 


251 (-) 



(a) Ref. H 

(b) Ref. [33 

Note that 5=1 and 5/2 are used for the fits of Mk(T) while Tq and T crossover (T cr ) are determined from the 
experimental resistivity curves. The T cr 's observed in Gao.94oMn .o6oAs (annealed: 370°C) and Gao.g3oMn .o7oAs (as 
grown) are 10 K and 12 K, respectively, which are very close to the calculated values of 8 K and 12 K, respectively 
from Eq. (fTOjl . The calculated carrier density is 2 x 10 20 cm -3 , using Ej =p Ep ~ 10 K, and Eq. ||7J). In this 
calculation, we take the effective mass as m* = 10mo, in accordance with the optical spectroscopy measurements [23| 
and compensation due to annealing, where mo is the rest mass of an electron. 

Figures HJc) and0|d) show the normalized magnetization, M a (T). M Pi td,k(T) are compared with the experimen- 
tally determined magnetization [8[ (M exp (T)) as depicted in Fig.[4|d). One can easily notice the inequality, Mtd{T) 
> M K {T) > M P (T) > A/ oxp (T) from Figs.g^c) anddjd). As such, M p (T)Js the best fit as compared with M cxp (T), 
better than the models developed by Tinbergen-Dekker [28| and Kasuya [29[ . 
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T(K) T(K) 

FIG. 4: (a) Calculated, using Eq. (|10[) , and experimental resistivity curves (p(T)) as a function of temperature for Gai-^Mn^As. 
BO and Al indicate Gao.93oMno.o7oAs (as grown) and Gao.94oMno.o6oAs (annealed at 370°C), respectively, (b) Calculated, 
using Eq. (|14|) . and experimental p(T) curves for annealed non-ferromagnetic Gao.930Mno.070 As samples. Bl and B2 indicate 
Gao.93oMno.o7oAs (annealed at 370°C) and Gao.g30Mno.070 As (annealed at 390°C), respectively, (c) and (d) show the T variation 
of calculated normalized magnetization M a (T) where a = K (Eq. l|13j l). TD (Eq. I|12|l ). p (Eq. exp (experimentally 

determined magnetization data) curves with spin quantum number S = 1 for x = 0.060 and 0.070 respectively. Mk(T) is also 
calculated with S = 5/2 as shown in (c) and (d). The experimental results M exp (T) for x = 0.070 (as grown) are shown in (d). 



B. Mn x Ge!- x 
1. p(T>T c ) 

Another DMS material that we will consider here is the Mn^Gei-^ system. It is a p-type DMS, with carrier 
density of the order of 10 19 — 10 20 cm~ 3 for 0.006 < x < 0.035 at room temperature [3J|. From the Ei model, 
Mn 3+ (18.910 eV atom" 1 ) substitution into Ge 4+ (25.941 eV atom" 1 ) sites will shift the p(T) curve downwards since 
{Ei) Ge 4+ > (£ , /) Mn 3+. This has been observed experimentally (see Fig. 2B of Ref. [34j]) where, p(T,x = 0.009) 
> p(T,x — 0.016) > p(T,x — 0.02). Subsequently, we can estimate the maximum valence state for the Mn ion in 
Mn^Gei.z using Eq. ^ and we find, S = 0.137 or Mn <3137+ : [^ / (Mn 4+ )] + i[£ / (Mn 3+ )+£ / (Mn 2+ )+£; / (Mn 1+ )] = 
i[£/(Gc 4+ ) + £;/(Ge 3+ ) + £j(Gc 2+ ) + Sj(Ge 1+ )]. Therefore, (51.200 x S) + 18.910 = 25.941 and 5 = 0.137. If the 
valence state of Mn is larger than 3.137+, then the p(T,x) curve is expected to shift upwards with x. The result, 
p(T, x = 0.02) « p(T, x = 0.033) indicates that Mn 3+ in this doping range may not have substituted Ge, instead it 
could have occupied non-substitutional sites that eventually causes the change in the valence state of Mn. Since the 
resistivity, p(T,x = 0.02) « p(T,x — 0.033), we can calculate the maximum increment of Mn 4+ content from x = 
0.02 to x = 0.033. Using S = 0.137, we estimate that there are 13.7% more Mn 4+ in Mno.033Geo.966 as compared to 
Mno.02Geo.98 samples. 
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2. p{T > 0) 



Resistivity measurements 34| and the fit using Eq. (fT0|) are shown in Fig.[5]Ja). From the fit, we find that Ej =F Ep 
— 15 K for Mno.02Geo.98 and the hole density is 2.38 x 10 19 cm" 3 using Eq. 0, where we take m* h = m , as the 
lower limit. This lower limit value is still comparable with the experimental value given above. We also find that 
the semiconductor-like behavior of p(T, x — 0.02) below Tc is not exponentially driven as indicated by Eq. (fTU|) in 
Fig. Eta). 




FIG. 5: (a) Calculated, using Eq. (|10p , and experimental resistivity curve p{T), for Mno.02Geo.9s- Also shown is the calculated 
p(T) curve with the additional constraint, dM p (T)/dT — in Eq. I|10p to emphasize the influence of M P (T) for an accurate 
fit. For these calculation, A — 25, B = 1060 and Ei + Ef = 15 K. We also show that the T-dependence of p(T) as calculated 
using Eq. (|14|l is in accordance with charge current (J e ) only. This result lacks the ability to capture the experimental result. 
In this case, we used A = 1.8 and Ei + E F = 15 K. Both Ei + E F = 15 K and Ei + E F = 80 K give p of the order of 10 19 cm" 3 
using Eq. Q and m* h = mo (rest mass), (b) shows the calculated and experimental T dependence of the magnetization curves 
M a (T) where a = K (Eq. (|13[) h P (Eq. (|ll[) h ex P (experimentally determined magnetization data). Notice the inequality, 
M P (T) > M exp (T) that arises as a result of the principle of least action. 



The pronounced effect of Eq. (fTTj) can be noticed by comparing the calculated plots between Eq. (fTO)) . and Eq. (fT0|) 
with the additional constraint dM p (T)/dT = 0, as depicted in Fig. EJa). The normalized magnetization, Mk. p . kp(T) 
for Mno.02Geo.98 is given in Fig.EJb). Again, we find that M p (T) gives the best fit for the experimental data, which is 
better than Mk(T) and Mtd(T). However, M P (T) is significantly larger than M oxp (T), which makes the fit poor. The 
reason is that the resistivity measures only the lowest Ej path, regardless of the temperature and with easily-aligned 
spin path that gives rise to high conductivity and complies with the principle of least action. That is, the ability of 
both J e and J se to follow the easiest path. On the contrary, magnetization measurements quantify the average spin 
of an ensemble of electrons, which will be smaller in magnitude as compared with transport measurements. 



C. Lai -xCajMnOs 



1. p(T>T c ) 

For manganites, the Lai-^Ca^MnOa system has a maximum Tc of 260 K (Ref. [35|). Unlike DMS, the resistivity of 
this class of materials above Tc is exponential and thus gives rise to the large drop in resistance below Tc that leads 
to the CMR effect. We will evaluate this scenario with the ionization energy concept and show the validity of iFDS in 
both DMS and manganites. Based on Ej model, Ca 2 + (£7 = 9.000 cV atom" 1 ) < La 3+ (£/ = 11.940 eV atom" 1 ), 
therefore the resistivity p{T, x) curve is expected to shift downward with Ca 2+ doping. On the contrary, the overall 
p(T) curve between x = 0.1 and 0.2, above Tc are almost identical (see Fig. 3 of Ref. |35||). Again, this could be due to 
the change in the valence state of Mn, as a result of Mn, Ca and/or La occupying the non-substitutional sites, provided 
that the valence states of Ca 2+ and La 3+ are invariant to doping and defects. Indeed, the content of Mn 4+ is found to 
increase with Ca doping [35|, from 19% (for x = 0.1) to 25% (for x = 0.2). We can calculate the maximum increment 
of Mn 4 + from x = 0.1 to 0.2, and we obtain, 5.75% using Eq. (TU): [5 x £;(Mn 4+ )] + i[£ 7 (Ca 2+ ) + £/(Ca 1+ )] = 
i[£:/(La 3+ ) + £;/(La 2+ ) + £/(La 1+ )]. Therefore, (51.200 x 5) + 9.000 = 11.940 and 8 = 0.0575. This value is 
remarkably close to the experimental value of 6%(=25%— 19%), determined via redox titrations [35l |. This implies 
that, E l j cal oc Ei, which in turn implies that the averaged many-body or lattice potential (j3) can be approximated as 
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a constant. This comes as no surprise because the averaged crystal or lattice potential is indeed a constant due to 
the periodicity. 

2. p(T > 0) 

Using Eq. (jTUJ) for La 1 _ I ,Ca a; Mn03 [35[, Ej =F Ep is calculated for x — 0.1 and 0.2 samples and we obtain 0.121 
eV (1400 K) and 0.112 eV (1300 K) respectively. The calculated carrier density, using m* = mo and Eq. {7J gives 
10 17 cm -3 . In the presence of the magnetic field, H = 6 Tesla, we obtain Ej =p Ep = 0.0776 eV for x = 0.2 and its 
hole concentration p = 10 18 cm -3 . The fits are shown in Figs.[6ja) and[U(b), while the fitting parameters are listed 
in Table [U The 5 value determined previously (5.75 %) needs to be corrected because in the previous calculation, we 
used (Ei =p Ef) x= o.i ~ (Ei =F Ef)x=o.2 that gives maximum 5. Therefore, the actual increment after the fitting is 
6 = 5.75 x (0.112eV/0.121eV) = 5.32%. 




150 150 300 



T(K) T(K) 

FIG. 6: (a) shows the experimental results for the resisitivity versus temperature, p(T), for Lai-^Ca^MnOa for x = 0.1, 0.2 
and 0.2 (6 Tesla) and (b) the calculated p{T) using Eq. (|10p . All fits are indicated with solid lines, (c) and (d) show the T 
variation of the magnetization, M a (T), where a — K (Eq. {T3}) and p (Eq. (|11|) ') with S = 1 for x = 0.1 and 0.2, respectively. 
The experimental M oxp (T) curve for x = 0.2 is given in (d). 

Figures [6^c) and(6^d) depict the calculated M a (T) with S = 1 and M exp (T) for x = 0.2, respectively. It can be 
seen that the magnetization curve calculated from Eq. (|11| shows better agreement with the experimental data than 
Mk(T). Hence, the ionization energy model is suitable for both types of ferromagnets, be it diluted or concentrated. 
However, this does not imply that the ferromagnetic interactions are identical between diluted and concentrated 
ferromagnets. 
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IV. CONCLUSIONS 



In conclusion, we have developed a theoretical model based on the ionization energy concept that can be used to 
analyze the evolution of the resistivity versus temperature curves of both diluted and concentrated ferromagncts, for 
different doping elements. By identifying the cause that deviate the ionization energy from its averaged value, we 
come to understand how defects and clustering contribute to the changes in valence states of ions and eventually how 
they affect the resistivity of ferromagnets. 
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V. APPENDIX 



A. Additional Notes 



The relation between the electron affinity and ionization energy is through r, where the electron affinity is connected 
to r — > oo, and for any other finite r, electron affinity is undefined microscopically and it is not a good quantum 
mechanical variable, as opposed to our £. For example, if we were to work with the electron affinity then we will face 
two problems, firstly, considering the anions, electrons from different cations will be treated with equal affinity toward 
the anion, which is incorrect. Secondly, considering the cations, there is no mathematical analogy as the first electron 
affinity, second, etc., which we have used for first, second, etc. ionization energies. Hence, the electron affinity in 
our case, is not a good many-body variable, except for qualitative descriptions. As for the reason behind the name, 
ionization energy is that in the early stages of the ionization energy theory, we used the atomic ionization energy 



(r — > oo) as the input parameter to compute carrier concentrations 17 1, and for this reason, £ was labeled as the 
ionization energy. The technically correct label is the excitation energy. We can also explain the resistivity above the 
Curie temperature with any exponential function that either contains the activation energy or the energy related to 
the electrons hopping rates. However, the whole point of our theory is to associate our exponential function (as given 
in Eq. (|10[) ) with each constituent atoms that exist in a given compound. This enable us to estimate and predict 
the evolution of the carrier density and the resistivity, before the experiments are carried out. On the contrary, the 
theories with activation energy and with hopping rates do not give us that freedom, hence these other approaches 
cannot predict the evolution of the resistivity for different doping a priori. In other words, these other approaches 
need to estimate the carrier density by other means to feed into their theory so as to predict the doping dependent 
resistivity. Apart from that, our resistivity theory captures the resistivity curve completely, from T > to T > Tq 
with a single equation, as presented by Eq. (jTUJ) . We do not have one equation for T > Tc, another for T < Tc 
and another one for T < T cro8Sove r- Our theory captures all three mechanisms with a single equation, which will be 
very useful for experimental evaluations. However, this approach has not been developed to the extend where we can 
apply it to magnetic heterostructures to estimate the spin transmission probability as discussed by Egues [36[ , Dai et 
al. [37|, and Papp and Peeters [38l. l39l. liol] . 

The proof that is given in the subsequent section is a form of proof-of-existence. This appendix is to prove that 
Eq. ((T|) is mathematically and physically correct, regardless of the type of potential used. In order to prove this, we 
used the simplest case of ID harmonic oscillator for convenience. Even if we were to use the hydrogenic wavefunction, 
we will still arrive at Eq. (fT]) as derived in Ref. [22| . In mathematics, there are many techniques to prove the existence 
of a solution for a differential equation. Of course, one of the techniques is to solve the equation. There are also the 
so-called existence and direct proofs. Now, to prove the validity of Eq. (fT]) is quite easy because £ is also an eigenvalue 
and we did not touch the Hamilton operator. Hence, we are not required to solve Eq. ([1]) in order to prove its validity. 
In other words, we can prove by means of constructive (existence) and/or direct proofs, by choosing a particular form 
of wavefunction with known solution and then calculate the total energy by comparison. In doing so, we will find 
that the total energy is always given by Eq ± £, as it should be. For example, for ID harmonic oscillator, the known 
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solution is a gaussian function, for ID Dirac-dclta potential, the known solution is an exponential function, for a ID 
square well potential, the known solution is a sinusoidal function, for a 3D Coulomb potential in hydrogen atom, the 
known radial solution is an exponential function. Therefore, it depends on which system we are interested in and if 
we are interested in ID harmonic oscillator, then we write the solution for Eq. ([T]) in its gaussian form and then we 
derive the Schrodinger equation and its wavefunction in terms of Eq ±£. After that, we compare this new exponential 
wavefunction with the known solution for ID harmonic oscillator. When we compare the constants, we will be able 
to show that the total energy is always given by Eq ± £. 



B. ID harmonic oscillator: ionization energy as the eigenvalue 

The one-dimensional Hamiltonian of mass m moving in the presence of potential, V(x) is given by (after making 
use of the linear momentum operator, p = —ihd 2 /dx 2 ) (4l| . 



2m dx 2 



ip = Eip, 

= (Skin + V pot )<p, (17) 



where E, £kin, V(x) and V po t denote the total energy, kinetic energy, potential energy operator and the potential 
energy, respectively. 
We define. 

± £ := Skin- E + V pot , (18) 

such that ±£ is the energy needed for a particle to overcome the bound state and the potential that surrounds it. 
The + sign for ±£ is for the electron (0 — > +oo) while the — sign is for the hole (— oo — > 0). i?kin and Eq denote the 
total energy at V po t = and the energy at T = 0, respectively, i.e., E^ n = kinetic energy. In physical terms, £ is 
defined as the ionization energy. That is, £ is the energy needed to excite a particular electron to a finite distance, r, 
not necessarily r — > co. 

On the other hand, using the above stated new definition (Eq. (fT5]l ) and the condition, T = and V po t = 0, we 
can rewrite the total energy as 

E = E kia = Eq ± e (19) 
From Eq. (fT5|) we have E = Eq ± £ = E kin + V pot , therefore 



H<p=(Eo±S)<p, 



(20) 



where the total energy is given by E = Eq ± £ and (p is the ionization energy based wavefunction while H is the 
Hamilton operator. 

Proof: Assume a solution for Eq. (|2"D|) at n = state (ground state) in the form of (p„-q(x) = Cexp[— ax 2 } in 
order to be compared with the standard harmonic oscillator wavefunction [4l| . 

(p (x) 





1/4 




raw 




irh 


exp 


_~ ~2h x 



(21) 



where H = h/2ir, h is Planck's constant and to is the frequency. Therefore, we obtain lu.ip(x) = InC — ax 2 hie, 
= —2ax and 9 = 2aip(x)[2ax 2 — 1]. On the other hand, we can rewrite Eq. (|2"0|) to obtain 



d 2 (f 2m 

— = -—[ Eo ±^, 



(22) 



where E and Eq for a given system range from +oo to for electrons and to — co for holes which explains the ± 
sign in £. Using Eqs. (fT7)l . (|2U)) and, (|2"2"j) . we obtain a = ^(Eq ± £). Normalizing ip gives 
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1 = 



+ 00 



\ip(x)\ 2 dx = 2\C\ 



~ 2ax dx 



2icr- 



7T 

2a 



1/2 



(23) 



and C 



2m{E ±j) 



1/4 



. Consequently, 



ip(x) = 



~2m(E ±£)' 


1/4 


exp 


irh 2 



(24) 



Now, we compare Eq. (|24|) with Eq. (|2~Tj) . In doing so, we can show that the ground state energy is, 

1 



either from equating Imw/nHl ^ 4 = [2m(E ± £))/7rfr 2 ] 1/4 or —mx 2 uj/2h = —m(Eo ± ^)x 2 /h 2 . 

C. Expectation value for V(x) 

In this section, we will show that the total energy, Eq ± £ is a function of the potential energy That is, using 
Eqs. (|17p . (fl~8|) and (f20|) . we will show the potential energy can be written in terms of Eq ± £. For example, from 
Eq. (|2~4"|) the harmonic oscillator Schrodinger Eq. can be shown as 



\hu = E ± e, 



1/4 



(25) 



h 2 d 2 v 

2m dx 2 

Therefore, the potential energy is given by 



2m , .to, 
- JiT (EQ±0 2 x 2 ^(EQ±0 



(26) 



V(x) = ^(Eq±H) 2 x 2 . 



(27) 



Proof: From Eq. (j2~7| , we can write 



H 



2m 



2m 

1 — (Eq ±£)x 



2 1 



(28) 



Therefore, the ladder operator can be written as 



a± 



= A 



. 2m , „ 
Tt p+ ( _(£ 0± £) ) x 



(29) 



t4 is a factor that will be used to derive the expectation value of V(x). Taking K = =j^(Eq ± £), we obtain 

a± = A[^if> + Kx}. 
Since the commutation relation [4lj . [x,p] = iK, then 



(30) 



Using Eq. pg]). we get 



[a-,a+] = 2A 2 Kh = -[a+,a_]. 



(31) 
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a_a + = A 2 (2mH + KK), 

1 ra T a+ 
H = - — -^p- =f • 



Consequently, we can show that 



2m 



1 



= 2^ 



A 2 



a+a_ 



(32) 
(33) 



A 2 

(E ± 0*>- 



(a+ip) 



(34) 



Subsequently, H(a_tp) = [E ±£ — 2(E ± £)] a_(y9. Applying the condition [4JJ for the ground state, </>„=o such 
that 



will lead us to 



a_v? Il=0 = 0, 



E n = ip n = — Hlp„ = o — (Aq ± £)<p rl= o 



(35) 



(36) 



Recall that Eq is the energy at T = and A n= o denotes the energy for the n = state. Finally, utilizing Eqs. |36|) 
and (|34[) . we obtain 



A„ = (l + 2n)(£o±0- 
Subsequently, we find that (using Eqs. ([20]), d33j) and IpTF])) 



(37) 



a + a-(p n = n2A 2 Khip n , 

a-a + Lp n = (n + 2A 2 AA)< y 5„ [o_, a+] = 2A 2 AA 



Now, using the identity [4l| 



We find 



(f*)(a±g)dx= (a T f*)(g)dx, 



(a + (p n )* (a + ip n )dx 



(a-a + ip n )*ip n dx 

poo 

(n + 2A 2 Kh) / |^„| 2 dx, 



(a_^n)*(a_y> n )dx 



(a + a_^„)*<^„(ia; 

) 

/OO 
|^„| 2 dx. 
-oo 



On the other hand, we have 



a±(f n (X <p n ±i a±ip n = 7„ if n ±l, 
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where 7+ and j n are the proportionality factors, which can be determined from 

/oo 
\tp n ±i\ 2 dx, 
-00 

hence, we can now write a+<p n = (n + 2A 2 Kh) 1 / 2 Lp n+ i and 

a-ip n = (n2A 2 Kh) 1 / 2 tp n -i. (38) 

We can rearrange Eq. (j30|) to get 

2 _ (a+ + a-) 2 f x 

As a consequence, (from Eqs. (|3T))) and (|3"9")) . also invoking the orthogonality in the last step) 

' 2m 

2A 2 Kh(l + n) + n 



(V(x)) = (^(E ±O 2 x 2 



8Ah 



From Eq. (|37|) . we know that 



(40) 



(V(x)) = ±(l + 2n)(E ±Q, (41) 
the other half is due to kinetic energy [4l|. Putting Eqs. (|4"0")) and (|4"Tj) together leaves us with 

2A 2 Kh(l + n) + n \, 

^-^ = -(^^±0, (42) 

Therefore A 2 = 27^- Hence, the commutation relation given in Eq. pip can be rewritten as 

[a-,a + ] = 2A 2 Kh = 1 = -[o+,a_]. (43) 
As a result of this, indeed the potential energy is given in terms of E ± £ from Eq. (|4Tj) . 



[1] I. Zutic, J. Fabian and S. D. Sarma, Rev. Mod. Phys. 76 (2004) 323. 

[2] H. Munekata, H. Ohno, S. von Molnar, A. Segmuller, L. L. Chang and L. Esaki, Phys. Rev. Lett. 63 (1989) 1849. 
[3] R. von Helmolt, J. Wecker, B. Holzapfel, L. Schultz and K. Samwer, Phys. Rev. Lett. 71 (1993) 2331. 
[4] M. Eginligil, G. Kim, Y. Yoon, J. P. Bird, H. Luo and B. D. McCombe, Physica E 40 (2008) 2104. 
[5] I. Kanazawa, Physica E 40 (2007) 277. 

[6] L. Saminadayar, P. Mohanty, R. A. Webb, P. Degiovanni and C. Bauerle, Physica E 40 (2007) 12. 
[7] T. Jayasekera, N. Goel, M. A. Morrison and K. Mullen, Physica E 34 (2006) 584. 

[8] A. van Esch, L. van Bockstal, J. de Boeck, G. Verbanck, A. S. van Steenbergen, P. J. Wellmann, B. Grietens, R. Bogaerts, 

F. Herlach and G. Borghs, Phys. Rev. B 56 (1997) 13103. 
[9] T. Omiya, F. Matsukura, T. Dietl, Y. Ohno, T. Sakon, M. Motokawa and H. Ohno, Physica E 7 (2000) 976. 
[10] J. M. Luttinger and W. Kohn, Phys. Rev. 97 (1955) 869. 

[11] T. Jungwirth, J. Sinova, J. Macek, J. Kucera, and A. H. MacDonald, Rev. Mod. Phys. 78 (2006) 809. 

[12] E. H. Hwang E H and S. D. Sarma, Phys. Rev. B 72 (2005) 35210. 

[13] M. P. Lopez-sancho and L. Brey, Phys. Rev. B 68 (2003) 113201. 

[14] C. Sen, G. Alvarez, H. Aliaga and E. Dagotto, Phys. Rev. B 73 (2006) 224441. 



17 



M. Mayr, A. Moreo, J. A. Verges, J. Arispe, A. Feiquin and E. Dagotto, Phys. Rev. Lett. 86 (2000) 135. 
T. Dietl, Physica E 35 (2006) 293. 

A. D. Arulsamy, Physica C 356 (2001) 62, |arXiv:cond-mat/0402153| 
A. D. Arulsamy, Phys. Lett. A 300 (2002) 691. 

A. D. Arulsamy in Superconductivity research at the leading edge (ed Lewis, P. S.) 45 (Nova Science Publishers, New York, 
2004). 

A. D. Arulsamy, Phys. Lett. A 334 (2005) 413. 

P. W. Anderson, Science 177 (1972) 393 . 

A. D. Arulsamy, |arXiv:physics/0702232^ 9; larXiv:0807.0745l 

K. S. Burch, D. B. Shrekenhamer, E. J. Singley, J. Stephens, B. L. Sheu, R. K. Kawakami, P. Schiffer, N. Samarth, D. D. 
Awschalom and D. N. Basov, Phys. Rev. Lett. 97 (2006) 87208. 

O. D. D. Couto, J. Rudolph, F. Likawa, R. Hey and P. V. Santos, Physica E 40 (2008) 1797. 

B. Kaestner, J. Wunderlich, T. Jungwirth, J. Sinova, K. Nomura and A. H. MacDonald, Physica E 34 (2006) 47. 

H. A. Bethe and E. E. Salpeter, Quantum mechanics of one- and two-electron atoms (Springer- Verlag, Berlin, 1957). 
S. Murakami, N. Nagaosa and S. C. Zhang, Science 301 (2003) 1348. 
Tineke Van Peski-Tinbergen and A. J. Dekker, Physica 29 (1963) 917. 
T. Kasuya, Prog. Theor. Phys. 16 (1956) 58. 

F. Matsukura, H. Ohno, A. Shen and Y. Suga wara, Phys. Rev. B 57 (1998) R2037. 
M. J. Winter <http://www.webelements.com>. 

A. Oiwa, S. Katsumoto, A. Endo, M. Hirasawa, Y. Iye, H. Ohno, F. Matsukura, A. Shen and Y. Sugawara, Solid State 
Commun. 103 (1997) 209. 

P. Mahadevan and A. Zunger, Phys. Rev. B 68 (2003) 75202. 

Y. D. Park, A. T. Hanbicki, S. C. Erwin, C. S. Hellberg, J. M. Sullivan, J. E. Matson, T. F. Ambrose, A. Wilson, G. 
Spanos and B. T. Jonker, Science 295 (2002) 651. 

R. Mahendiran, S. K. Tiwary, A. K. Raychaudhuri, T. V. Ramakrishnan, R. Mahesh, N. Rangavittal and C. N. R. Rao, 

Phys. Rev. B 53 (1996) 3348. 

J. C. Egues, Phys. Rev. Lett. 80 (1998) 4578. 

N. Dai, H. Luo, F. C. Zhang, N. Samarth, M. Dobrowolska and J. K. Furdyna, Phys. Rev. Lett. 67 (1991) 3824. 

G. Papp and F. M. Peeters, Phys. Stat. Sol. (b) 241 (2004) 222. 
G. Papp and F. M. Peeters, Appl. Phys. Lett. 78 (2001) 2184. 
G. Papp and F. M. Peeters, Appl. Phys. Lett. 79 (2001) 3198. 

D. J. Griffiths Introduction to quantum mechanics (Prentice-Hall, New Jersey, 1995). 



